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Abstract 


Regular structures generated by random interactions in energy cen¬ 
troids defined over irreducible representations (irreps) of some of the 
group symmetries of the interacting boson models sdIBM, sd( 7 lBM, 
sdlBM-T and sdIBM-ST are studied by deriving trace propagations 
equations for the centroids. It is found that, with random interactions, 
the lowest and highest group irreps in general carry most of the prob¬ 
ability for the corresponding centroids to be lowest in energy. This 
generalizes the result known earlier, via numerical diagonalization, for 
the more complicated fixed spin (J) centroids where simple trace prop¬ 
agation is not possible. 
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Johnson, Bertsch and Dean in 1998 [1], using the nuclear shell model found that 
random two-body interactions lead to ground states, for even-even nuclei, having 
spin O’*" with very high probability. Similarly, Bijker and Frank [2] using the inter¬ 
acting boson model with s and d bosons (sdIBM) showed that random interactions 
generate vibrational and rotational structures with high probability. Further studies 
using the shell model, fermions in one or two j orbits, sd, sp and sdg IBM’s, bosons 
in a single i orbit etc., revealed statistical predominance of odd-even staggering in 
binding energies. O’*", 2~^, 4+, ... yrast sequence, regularities in ground states in par¬ 
ity distributions, occupation numbers and so on; see [3-10] and references therein. 
Notably, Zelevinsky et al [4] introduced the idea of geometric chaos as a basis for 
the regularities observed in shell model studies. Similarly, Zhao et al [8] developed 
a prescription based on sampling of the corners of the parameter space and Bijker 
and Frank [7,9] employed mean-held methods. The unexpected results for regulari¬ 
ties with random interactions are reviewed in [6,10]. As Zhao et al stated [10]: ”a 
more fundamental understanding of the robustness of 0^^ dominance is still out of 
reach”. Therefore, going beyond the ground states and near yrast levels, energy cen¬ 
troids, spectral widths and correlations among them are also being investigated by 
several groups [4,5,11-13] as they are expected to give new insights into regularities 
generated by random interactions. For example, Zhao et al [11,12] initiated the study 
of energy centroids and analyzed £xed-L (£xed-J, JT) centroids in IBM’s (in shell 
model) spaces. They found that (or Jmin) and Lmax (or Jmax) will be lowest 

with largest probabilities and others appear with negligible probability. Similarly Pa- 
penbrock and Weidenmiiller [13] recently analyzed the structure of hxed-J spectral 
widths for fermions in a single-j shell. 

An interesting and important question is the extension of the spin zero ground 
state dominance (and also other regular structures seen in shell model and IBM 


2 



studies) to group theoretical models with hamiltonians preserving a symmetry higher 
than J (or L). Similarly one may consider centroids and variances dehned over 
good or broken symmetry snbspaces. They open a new window to the regularities 
of many-body systems in the presence of random forces. Initiating work in this 
direction [9], recently random one plus two-body hamiltonians invariant with respect 
to 0(A/i) © 0(A/2) symmetry of a variety of interacting boson models are used to 
investigate the probability of occurrence of a given (ci;ici; 2 ) irreducible representation 
(irrep) to be the ground state in even-even nnclei; [cui] and [uj‘ 2 \ are symmetric irreps 
of 0(A/i) and 0(A/2) respectively. It is fonnd that the O"*" dominance observed in 
ground states of even-even nuclei extends to gronp irreps. The pnrpose of this paper 
and others to follow is to go beyond this and stndy regnlarities, generated by random 
interactions, in energy centroids, variances etc. dehned over group irreps. Reported 
in this rapid commnnication are the results of a hrst analysis of energy centroids with 
examples from sdIBM, sdgflBM, sdIBM-T with the bosons carrying isospin (T) and 
sdlBM-ST with the bosons carrying spin-isospin {ST) degrees of freedom. Before 
proceeding further, it is important to stress that energy centroids (also variances) 
can be calculated as a function of particle number (m) and the qnantum numbers 
labeling the group irreps, without recourse to the construction of the hamiltonian 
matrix. The principle used here is trace propagation, a snbject introdnced in the 
context of statistical nuclear spectroscopy by French [14,15]. Readers not interested 
in the details of group algebra and derivation of trace propagation eqnations for the 
energy centroids (given by Eqs. (5), (6), (7) and (8)), may leap ahead to the discussion 
of results starting just after Eq. (8). 

Let us begin with the spectrum generating algebra (SGA), say Gi, of a group 
theoretical model with all the many particle states in the model belong to the irrep 
Ti of Gi- For example, the SGA Gi for sdIBM is U{6). Now the average of an 


3 



operator 0{k) of maximum body rank k over the irreps r 2 of a subalgebra G 2 of Gi 
{G 2 in general denotes a set of subalgebras contained in Gi and r 2 denotes all their 
irreps) is dehned by 


^ ^ {T^(3T2a I 1 I T^(3T2a) 

/9 a&V2 

( 1 ) 

In Eq. (1), (3 labels the multiple occurrence (multiplicity) of r 2 in a given Ei irrep (i.e. 
in the reduction of Ei to E 2 ). Removing the denominator in Eq. (1) gives the trace 
over (Ei,E 2 ) space, i.e. tr[0{k)Y^’^^■ General theory for propagation of traces of 
operators over irreps of group symmetries is developed in Refs. [16-19]. In particular, 
Quesne [16] showed that, for Gi D G 2 , trace propagation over the irreps Ei and E 2 of 
Gi and G 2 algebras is related to the so-called integrity basis of G 2 in Gi which gives 
the minimal set of G 2 scalars in Gi. As discussed in Refs. [17,18], it is seen that in 
general the multiplicity of E 2 in a given Ei irrep results in the propagation of matrix 
of traces tr[0{kyf^^^'^^' = Ha (ri/5E2Q; | 0{k) \ Tij3'T2a). However, quite often the 
trace of this trace matrix or its average, as given by Eq. (1), which is important in 
applications, may not propagate in a simple manner. There are approximate methods 
for propagating trace of the trace matrix and they are signihcant in particular when 
the integrity basis contains far too many operators [18]. A very important example 
here is £xed-L averages in IBM’s or fixed-J (and JT) averages in the shell model. For 
these, it is not possible to write a simple propagation equation in terms of the defining 
space averages. On the other hand traces over irreps of group symmetries (higher than 
J symmetry) can be propagated in many situations using Casimir invariants. French 
and Draayer [19] showed that by simple counting of irreps of G 2 in Gi and the scalars 
one can construct in terms of the Casimir invariants of Gi and G 2 will immediately 
confirm if propagation via Casimir invariants is possible; in this situation the integrity 


{C>(i.))'’..r> = y (ri^jr^a | o{k) \ ri/jr^o) / 

(3 a&T2 
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basis reduces to Casimir operators of Gi and G 2 - In this paper we restrict ourselves 
to examples in IBM’s where this result applies; Refs. [20,21] give first IBM examples. 

For IBM’s the SGA, called Gi above, is with A/” = 6 for sdIBM, 15 for 

sdgIBM etc. and its irreps Fi are labeled uniquely by the boson number m as all 
m boson states are symmetric with respect to U{M). Now, consider the average of 
an operator 0{k) over the irreps (m, F 2 ) with F 2 ’s being the irreps of a subalgebra 
G 2 of U{Af). Say the the number of (m, F 2 )’s, called F*’s hereafter, foTm<k is 
r. Also assume that there are r number of invariants Ci, i = 1,2,... ,r of maximum 
body rank k constructed out of the products of m and the Casimir invariants of G 2 . 
Then, for any irrep F°, clearly {O)^ = where a* are constants. The ads 

can be determined by assuming that the averages {O)^^ are known for the irreps Fb 
j = 1, 2,..., r. For example, F^’s can be chosen to be the irreps (m, F 2 )’s for m < k. 
With this, dehning the row matrices [C] and [Omp] and the r x r matrix [X] as 

[C] = , [0,np] ^ 0,np. = {Of , [X] X,, = , (2) 

the propagation equation is 

{Of = |C] |X|-' [oZ,] (3) 

As the eigenvalues of the Casimir invariants of the algebras U{M), 0{M) etc. are 
known, construction of [C] and [X] is easy. In the reminder of this paper the H 
is assumed to be (1 + 2)-body. As an example let us consider SU{3) centroids in 
sdIBM. Here Gi = U{6) and G 2 = SU{3). Simple counting of scalar in terms of 
the number operator n and the quadratic Casimir operator f and the cubic Casimir 
operator C 3 of SU(3), conhrm that they exhaust all the scalars needed for propagating 
{0{k)f’^^^^ for any k [18,21]. Note that (A/i)’s denote SU{3) irreps. Propagation 
equation for the energy centroids over SU{3) irreps can be written as _ 

ao + aim + a 2 m? + a^C 2 {\ij) where 
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( 4 ) 


C2(Ay(i) — I C 2 I (A/i)®^ — A^ -|- -|- AyU -|- 3(A -I- yu) 

Using Eqs. (3) and (4), the propagation equation, in terms of the energy centroids 
for m < 2, is [21] 

= i(2-3m + m2)(if)°Aoo)^(2m-m2)(hf)^420) 

+ [-|m + + iC 2 (Ay(i)] (hf) 2440 ) ( 5 ) 

+ \m + |m^ - ^C 2 (Ay(i) (if)^402) 

Eq. (5) extends easily to the SU{3) limit of p/IBM with U(10) SGA but not to 
sdg, sdgpf, etc. IBM’s. Now we will derive 3 new propagation equations for energy 
centroids. 

In the U (A/”) D J2i [U (M) ^ 0{N'i)] © symmetry limits of IBM’s, with the bosons 
carrying angular momenta £ 1 , £ 2 , • • • so that Ni = ( 2 U + 1 ) and Af = Afi, for a given 
ith U{Afi) generates number of particles m* in the orbit and 0{Afi) generates 

the corresponding seniority quantum number cUj. The number operators fii of U{Afi) 
and the quadratic Casimir operators of 0{Afi) or the corresponding pairing operators 
p2{0{Afi)) suffice to give hxed mu = {miUi,m 2 U 2 , ■ ■ ■) averages of H. Appendix 
A in Ref. [22] gives the explicit form of P2{0{Afi)) for a general situation. Fixed- 
mu centroids of H can be written as op' mpmj — Sij) + 

Hi Ci (P2{0{Afi))) . Solving for ai/s and q’s in terms of the centroids for m <2, 

the final propagation equation, for IBM’s with no internal degrees of freedom, is 

{H) = ^ miSi + ^ Vij mimj + - - - (U) ‘ 

i i>j i 

/T^\mi=2,i.Ji=0 _ iy\mi=2,ijJi=2 

+ - HM- - ~ o;j) (m* + cUj + A/i — 2) ; 

xy= mu, + p)] /(I + %)}-‘ E (2i +1) . 

L 

(U)”^*=2'-»=0 = ((UU)L. = 0 I U I = 0) , 

(^\^'rni=2,u)i=2 _ A/i(A/i+l) j A/i(Al+l) _ 
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Note that in Eqs. ( 6 ), are energies of the single particle levels with angu¬ 
lar momentum ii and = {{iiij)L \ V \ {ii£j)L) are two particle matrix el¬ 

ements of the two-body part of H. Also in Eq. ( 6 ), for s orbit rus = 2 and 

= 2 and there will be no two-boson state with Ug = 0. Eq. ( 6 ) for 
is given first in [23], i.e. for averages over the irreps of the algebras in the chain 
Usdg{l 5 ) D © [Ud{ 5 ) D Od(5)] © [Ug{ 9 ) D Og(9)]. Similarly Eq. ( 6 ) gives H 

averages over the irreps of UgdiQ) © Ud{5) D 0^(5) of sdIBM, Usdpfi^Q) © [Udi^) © 
Od(5)] © [Up{ 3 ) D Op(3)] © [Uf{ 7 ) D 0/(7)] of sdpfIBM etc. Moreover this extends 
easily (this will be discussed elsewhere) to IBM’s with internal degrees of freedom. Let 
us add that it is also possible to write down propagation equations for the variances 
{[H — results in [20,24]. 

In IBM-T with U{3J\f) D © [SUt{3) © Ot( 3)] where U{J\f) gives the spatial 
part (for sd, sdg, sdpf etc.) and 0^(3) generating isospin [25], it is possible to 
propagate the centroids ^ ^ Note that the U{J\f) irreps are 

labeled by {/} = {/i,/ 2 ,/ 3 } where /i > /2 > /s > 0 and m = /i + /s + /s- The 
corresponding SUt{3) irreps are (A,/i) = (/i — / 2,/2 — /s)- The SUt{3) to 0^(3) 
reductions follow from the formulas given by Elliott [26,27]. The scalars 1, h, h^, 
C2{SUt{3)) and and the energy centroids for m < 2, via Eqs. (2), (3) and (4), 
give 


(H) 




|m+^ 


(i7)0’(oo)-o + [2m - m2] 


-Am + 4^2 + lC2(A/i) - iT(T + 1)1 (H) 


2,(20),0 


+ 
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-|m + + ir(T + 1)] 

|lm+X-iC2(A,^)] . 

For sdIBM-T, starting with the general hamiltonian given in Appendix-A of [27] which 
contains the s and d boson energies and 17 two-particle matrix elements it is 

easy to write down m < 2] for m = 2 the two-boson isospins t uniquely 


(7) 
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define the corresponding SUt{S) irreps. Thus Eq. (7) for jg g^^gy apply 

for any m. 

In IBM-S'T, a group chain of interest is [28] U{6J\f) D U{J\f)<^[SUsT{(^) ^ Osr(6)] 
with U{Af) generating the spatial part and SUst{Q) [or Ust{Q)] generating spin- 
isospin part; note that the Wigner’s spin-isospin super-multiplet algebra SUst{4:) 
is isomorphic to Just as before, it is possible to propagate the centroids 

Here {/}’s are the irreps of U{J\f) or equivalently Ust{Q) and {/} = 
{/i,/ 2 ,/ 3 ,/ 4 ,/ 5 ,/ 6 } where fi = m and ft > fi+i > 0. The Osr(6) irreps are 
labeled by [a] = [cti, (J 2 , as] and the {/} to [a] reductions, needed for the results 
discussed ahead, follow from the analytical formulas given in [27] and the tabulations 
in [29]. Eqs. (2) and (3) give, using the quadratic Casimir invariants € 2^8 of Ust{^) 
and Osr(6), 


{H) 






+ 

+ 


^H-)OJoh[o] ^ [2m - m2] 

-|m + im 2 + JC2({/}) + i^C2([a])] (i7)2J2}.[2] 
+ i^C2({/}) - ^C2([a])] (i7)2J2}lo] 

fm +V-iC2({/})l ; 


( 8 ) 


where C 2 ({/}) = (C 2 (fJsT( 6 )))^'^^ = ^-f 7 - 2i) and C 2 ([a]) = 

i=l 

r 1 3 

(C 2 ( 05 t( 6 ))) = ^ a* (a* -h 6 - 2i). Diagonalizing C 2 {Ost{Q)) in the 1 ( 7172 )^ 5 ^) 

i=l 

basis and applying the resulting unitary transformation to the 77 matrix in this basis 
will give the input averages in Eq. ( 8 ). 

Now we will apply Eqs. (5)-(8) to study regularities generated by random interac¬ 
tions in energy centroids. In all the calculations used are independent Gaussian ran¬ 
dom variables with zero center and unit variance and a 1000 member ensemble. We be¬ 
gin with the simplest example of sdIBM centroids. The highest SU (3) irrep for a given 
m is (2m, 0) and Eq. (5) gives, ( 77 )'"^'^^^ — = [C 2 (A/i) — C2(2m, 0)] A/18 
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with A = Therefore the probability of hnding A to be positive or 

negative will simply give the probability for hnding the highest or lowest m particle 
SU{3) irrep to be lowest in energy. With the two-particle matrix elements chosen to 
be Gaussian variables (with zero center and unit variance), A itself will be a Gaussian 
variable with zero center. For m = 3k, 3k + 1 and 3k + 2, k being a positive integer, 
the lowest SU{3) irreps are (00), (20) and (02) respectively. They will be lowest in 
energy with 50% and the (2m, 0) irrep will be lowest in energy with 50% probability. 
Thus, it is easy to understand the regularities in centroids dehned over hxed SU (3) 
irreps in sdIBM with one plus two-body hamiltonians, without constructing the many 
boson hamiltonian matrix but just by using the propagation equation (5). 

In sdgfIBM, regularities in £xed-(ms, rud, Vd, nig, Vg) centroids are studied using the 
propagation Eq. (6). Ghoosing the 3 single particle energies {es,ed,eg) and the 16 
diagonal two-particle matrix elements with ii = 0,2 and 4 to be Gaussian 

variables, the probability for the centroid of a given (m^, m^, Vd, nig, Vg) configuration 
to be lowest is calculated for m = 6 — 25 and the results are shown in Fig. 1 
for m = 15. To maintain proper scaling, the e’s are divided by m and the by 
m(m — 1) just as in [2]. For the discussion of the results we dehne ti{x) such that 
■n{x) = 0 for a; even and 7r(a;) = 1 for a; odd. It is seen from Fig. 1, and also 
valid for any m, that the conhgurations {nis,nid = Vd = m — nig,nig = Vg = 0), 
{nis, md = ni- nig, Vd = 7r{nid),nig = Vg = 0), {mg, nid = Vd = 0,nig = Vg = m - nig) 
and {nis,nid = Vd = 0,nig = m — nig,Vg = n{nig)) exhaust about 91% probability. 
Moreover, the configurations with rrig = nid = 0 carry ~ 20%, nig = nig = 0 carry 
~ 21%, uig = m carries ~ 24% and m* % 0 but nid = 0 or nig = 0 carry ~ 26% 
probability. Thus the m* = m configuration and the four configurations with nig = 0 
are most probable to be lowest in energy. However, other conhgurations with nig % 
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0,m (they are 49 out of 1195, configurations in the m = 15 example) give non- 
negligible probability for being lowest. Thus about ~ 4% of the (m^, 
conhgurations will have probability to be lowest with random interactions. 

For sdlBM-T, it is easily seen from Eq. (7) that the one-body part of H will 
not play any role in the study of fixed-(A,/i)T centroids. Choosing to be 

Gaussian variables, the centroids are generated, using Eq. (7), for m = 10 — 25 
and for all allowed (A, fi)T. Some typical resnlts for the regularities are shown in 
Fig. 2. Firstly, for a given m the highest SUt{S) irrep is (m, 0) with Tmax = rn and 
Tmin = vr(m). For m = S/c, 3A;-|-1 and 3A;-|-2, with k being a positive integer, the lowest 
SUt{S) irreps are (00), (10) and (01) with T = 0, 1 and 1 respectively; for the later 
two situations the next lowest irreps are (02) and (20) respectively with = 0. 
For m = 3k, it is seen from Fig. 2 that the lowest SUt{3) irrep’s centroid (here T 
is nniqne) is lowest with ~ 35% probability. Similarly the highest irreps centroid is 
lowest with ~ 60% probability and this splits into ~ 30% each for the lowest and 
highest T’s. For m = 3k + 1 and 3k + 2, the probability for the centroid of the 
highest irrep to be lowest in energy is same as for m = 3k. However for the centroid 
of the lowest irrep, the probability is ~ 29% and the next lowest irrep appears with 
~ 6%. Thus in general the centroids of the highest and the lowest (for m = 3k + 1 
and 3k -|- 2, the lowest two) SUt{3) irreps exhaust about 95% of the probability 
for being lowest in energy. As the two particle centroids X* = qj-q linear 

combinations of E’s, it can be seen that they themselves are Gaussian variables. Note 
that _ [C 2 (A/i) - C 2 (m, 0)] Ai + [T(T + 1) - m(m + 1)] A 2 

where Ai = |A° -|- — |X^ and A 2 = |(X^ — X'^). Calculations with X*’s taken 

as Gaussian variables with same variance (actnally the variance of X° and X^ are 
same and that of X^ is ~ 20% higher) are carried out and it is seen that they give 
almost same results as in Fig. 2. 
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are de- 


For sdlBM-F'T, as seen from Eq. ( 8 ), the energy centroids 
termined by the 2 -particle averages ' and they 

are linear combinations of the two particle matrix elements in the \{iii 2 )LST) 

basis. Instead of choosing to be Gaussian variables, we have chosen, using the 

result found in the sdIBM-T examples, the three 2-particle averages to be Gaussian 
variables. Using this, the probabilities are calculated for various m values and some 
of the results are shown in Fig. 3. Firstly for a given m the highest {/} is {m}. The 
corresponding highest and lowest [a] are [m] and [ 7 r(m)]. For all m the centroid of 
the highest Ust{Q) irrep is lowest with ~ 56% probability and this splits into ~ 34% 
and ~ 22% for the highest and lowest OstG) irreps. For m = 6k, 6 /c ± 1, 6 A; ± 2 
and 6k + 3, with k a positive integer, the lowest Ust{6) irreps are those that can be 
reduced to the irreps {0}, {!}, ({2}, {1^}) and ({1^}, {21}) respectively. These irreps 
with the corresponding lowest [a] are lowest with probability ~ 43%. 

In conclusion, with random interactions, the lowest and highest group irreps (i.e. 
irreps of G 2 in Gi D G 2 ) carry most of the probability for the corresponding centroids 
to be lowest in energy. With the inclusion of a subalgebra (Gi D G 2 D G 3 ), these 
probabilities split into the probabilities for the corresponding lowest and highest irreps 
of the subalgebra. This is indeed the situation for all the examples discussed in this 
paper. Gontinuing with the process of embedding subalgebras, the 0(3) algebra 
generating L can be reached (with generalization for systems with LT, LST or JT). 
Then clearly the energy centroids of highest and lowest L’s should be most probable 
and this is found to be true numerically in [ 11 , 12 ]. An important aspect of the 
energy centroids is that they propagate via Gasimir invariants in many situations. 
New propagation equations are derived in this paper (Eqs. ( 6 ), (7) and ( 8 )). In fact 
there are many other situations where such equations can be derived; an example is 
for the centroids over the irreps [m^d(Asd/isd);^ 2 rp/(Ap//ip/)] of [Usd{6) D SUsdi^)] © 
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[Upf{10) D SUpf{3)] algebra of sdpflBM [30]. These will be discussed in a longer 
paper along with extensions of the present work to spectral variances and also to 
shell model symmetries. Finally, an important observation is that the propagators 
carry information about Gi D G 2 geometry (i.e. Gi D G 2 reduced Wigner coefficients 
and G 2 Racah coefficients) and thus it is plausible that propagation eqnations may 
be usefnl in qnantifying geometric chaos. This is being investigated and it should be 
remarked that only recently the role of Wigner-Racah algebra in two-body random 
matrix ensembles is established [31]. 

Thanks are due to Y.M. Zhao for useful correspondence and for making available 
Ref. [12] before it is snbmitted for pnblication. 
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FIGURES 



FIG. 1. Probabilities for sdgrlBM fixed-(ms, nid-, Vd, mg, Vg) centroid energies to be lowest 
in energy vs for a system of 15 bosons (m = 15). For each mg, the probability shown 
is the sum of the probabilities for the irreps with the seniority quantum number lowest 
(uf = 7r(m£)) and highest (ve = mg). Filled circles and stars are for configurations with 
md = 0 and mg = 0 respectively; they are joined by lines to guide the eye. Note that for 
rus = 15 both = 0 and mg = 0. 
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FIG. 2. Probabilities for the sdIBM-T’s (A/u)T centroid energies to be lowest in energy 
vs C 2 {\p)/rn?‘ for boson systems with m = 12,15,20 and 25. Except for the highest (A^), 
for all other (A/r)’s shown in the figure, the probabilities are for Tmax if A / 0 and /U 7 ^ 0 and 
they are for T^m if A = 0 or // = 0. For the irreps not shown in the figure, the probability 
is < 0.1%. All the points for a given m are joined by lines to guide the eye. 
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FIG. 3. Probabilities for the sdIBM-ST’s {/}[(t] centroid energies to be lowest in energy 
vs C 2 {{f})/rn? for boson systems with m = 6,7,9 and 10. For the structure of the irreps 
with probability > 2%, see text. For m = 7, 9 and 10 there is one additional irrep with 
~ 0.3% probability. For the irreps not shown in the figure, the probability is < 0.1%. All 
the points for a given m are joined by lines to guide the eye. 
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